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Nonlinear Transport and Current Fluctuation in an AB Ring with a Quantum Dot 

Junko Takahashi* and Shuichi Tasaki 
Department of Applied Physics, School of Science and Engineering, Waseda University, Tokyo 169-8555 

Nonequilibrium steady states are explicitly constructed for a noninteracting electron model 
of an Aharonov-Bohm (AB) ring with a quantum dot (QD) with the aid of asymptotic fields. 
The Fano line shapes and AB oscillations are shown to strongly depend on the bias voltage. 
Current fluctuations are studied as well. 
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1. Introduction 

For an electron transport in an Aharonov-Bohm (AB) 
ring with a quantum dot (QD), a coupling between the 
discrete states and continuum causes an asymmetry of 
the resonance peak in the conductance, 1 the so-called 
Fano resonance. 2 As well known, when one varies a mag- 
netic field piercing into the AB ring, the conductance 
oscillates (AB oscillation) . And the phase of the AB os- 
cillation changes by ir from one side of the Fano reso- 
nance to the other. 3 It is then interesting to see how 
these features appear in a far-from-equilibrium regime. 

Recently, the method of C*-algebra has been widely 
applied to rigorous investigations of nonequilibrium 
steady states (NESS) of infinitely extended quantum 
systems. 4-9 The method was originally developed to 
deal with equilibrium statistical mechanics of such sys- 
tems. 10 ' 11 Frohlich, Merkli and Ueltschi 9 have shown the 
existence of NESS and its characterization for junction 
systems, which consist of several free-fcrmionic reservoirs 
mutually interacting via bounded interactions. Their re- 
sults cover the AB ring with a QD in case that the elec- 
trodes are free fermion systems, but explicit expressions 
of observables at NESS are not given. 

In this article, exact NESS is constructed explicitly 
for a noninteracting-electron model of the AB ring with 
a QD. We then examine the bias voltage dependence 
of the Fano resonance shapes and the AB oscillation. 
Near the Fano peak, AB oscillation patterns are shown 
to be strongly affected by the bias voltage. This feature 
is shown to imply nonlinear conduction which is sensitive 
to the magnetic field. Although we consider a noninter- 
acting system, we believe that the Coulomb repulsion 
does not affect the essential features of the bias- voltage- 
dependence of the AB oscillation. The current fluctuation 
is investigated as well and we find that the interference 
effect appears in the fluctuation even when it cannot be 
seen in the average current. 

2. System and Nonequilibrium Steady State 

The system consists of two-dimensional electrodes and 
a quantum dot interacting with each other, which is 
schematically shown in Fig. 1. 12 Only the dynamics of 
electrons is considered and the electron-electron inter- 
action is neglected. Then, its Hamiltonian is given by 

*E-mail address: junchi@ruri.waseda.jp 



H = Hl + Hr + Hrj + Ht + Hlr , 

where H x = J dku).\a\ aX ak<TX (A = L (R)) is the en- 
ergy of the left (right) electrode and a ka x (A = L (Rj) 
is the annihilation operator of an electron with wave 
number k and spin a at the corresponding electrode. 
Throughout this article, the repeated spin index means 
the summation over it. The bias voltage V is taken into 
account as the difference of the one-particle energies: 
ujkn, = ^kL — eV\ The second term Hrj — £gc\c a stands 
for the energy of the QD, where e g is the dot level and 
Co- is the annihilation operator of the electron with spin 
a at the dot. The electron transfer through the QD is 
described by a tunneling Hamiltonian 

H T = J dk{u kL al aL c a + u kR al aR c a + (H.c.)> 
and the direct transfer by 

H L r = JJ dkdq{W kq e^al aL a qaR + (H.c.)}, 
where W kq — WukLUqR is assumed. The phase ip is pro- 
portional to the magnetic flux <j> piercing between the 
two paths: ip = eqj/ (he) with c the velocity of light. 
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Fig. 1. Schematic view of an AB ring with a quantum dot 

NESS can be simply characterized by in-coming fields 
Pka\- They are solutions of [H,pk<r\] = -^kxPkaX with 
the boundary condition: a kcr \(t) exp(iujk\t/h) — ► (3kaX as 
t —y — oo. When H does not have bound states, the orig- 
inal operators can be written as 



O-kaX = flkcrX + / dk! 



UkXUk'X 



Xk'xP, 



k'aX 



^k'X — ^>kX 
UkXUk'X 



iO A(w fc / A 

K kxPk'<rX 



dk 



u k >x ~ ^kx + iO H^k'x) 

UkL^kL^kaL , UkR^lkuPkaR 



(1) 



A(^fci) 
M\ (uJkx ) [W 2 (oJkx 



A(w/c_r) 

e g ) + 2Wcos<p] 



(2) 



where Xk\ = 1 

K feA = 1 + We^WfcA - € B ), n kX = l + We^M x {u kX ), 
A(x) = u(x)(x-e g )-J2\Mx{x)-2WcostpM L (x)M R (x) 
u(x) = 1 - W 2 M L (x)M R (x), and 
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M x (x) = / dk- 

J _ x ~ UJ kX_+ i0 
In the above, L = R, R = L, ipp — —<P and tpn = ip. 

NESS is constructed as follows: (1) Initially, the tun- 
neling interactions are turned off and the two electrodes 
are prepared to be in equilibrium with different tempera- 
tures (Tl and Tr) and chemical potentials (pp and fip). 
The initial dot state can be taken arbitrarily. Let po be 
the density matrix describing this initial state. (2) The 
tunneling interactions are turned on. Then, the whole 



state pt 



-iHt/h 



iHt/h b e g ms to evolve towards 



a steady state p+. One can actually show that, when 
t — > +oo, p t approaches to p + in a weak sense, namely, 
lim t ^ +00 ti(Apt) = tv(Ap + ) holds for every local observ- 
able (for detail, see Appendix). 

The so-obtained NESS p + satisfies Wick's theorem and 
is characterized by the two-point function: 

<#La/Wv> = F x (e)5(k - k')5 aa ,5 xy (3) 

where (...) = tr(...p + ) stands for the NESS average and 
F\(e) denotes the Fermi distribution function with tem- 
perature T x and chemical potential p x (A = L or R). 

Here, a remark is in order. In spite of the fact that the 
system is conservative (i.e., the evolution is governed by 
the Hamiltonian) , the state evolution exhibits 'dissipa- 
tion' in the weak sense. This behavior is due to the infi- 
nite extension of the system. Indeed, information about 
the initial state is 'dissipated' towards infinity (namely 
the Poincare recurrence cycle is infinite) and cannot be 
retrieved if we consider the average values of only the 
local observables. Thus, for every local observable A, 
the long-term limit lim t ^ +oc tr(Ap t ) exists. Note that 
the Keldysh Green's function method shares the same 
property when applied to the description of NESS. In 
reality, the electrodes are finite, but can be regarded 
as infinite systems compared with the mesoscopic sub- 
systems (AB-ring and QD). And the experimentally ob- 
served noncquilibrium states are considered to be well- 
described by NESS such as p + . 

Now we go back to the steady-state current. The op- 
erator J corresponding to the current from the left elec- 
trode is J = —jj- t J dk 



(J) 



if 



dku kL J\ 



dot 
k 



where j£ ot 
and J£™ = ie l 



kaL akaL and its average is 
+ e - J J dkdqWu kL u qR J^ (4) 

i{a\. aL c a ) + (c.c.) is the current via the QD 
(c.c.) is the current via the 



direct tunneling. With the aid of eqs.(l), (2) and (3), the 
current flowing from the left electrode is given by 



(J) = 



h 



deT(e){F L (e)-F R (e)} 



(5) 



where \i c — max(w L, ^or) an d T(e) is the transmission 
probability: 

Ak 2 u a D 2 



T(e) 



\A(e 



-{(W(e - e g ) + cos ip) 1 + sin 2 ip}. 



The constant D denotes the density of states of the 
two-dimensional electrodes. For the sake of the simplic- 
ity, the tunneling matrix elements u k L,UkR are taken 
to be symmetric and independent of the wave number: 



u k L = u kR = u. We remind that the bias voltage is equal 
to the chemical potential difference: pl — = e,V . In 
the following sections, we focus on the noncquilibrium 
properties at zero temperature: Tl — T R = 0. 

3. Transport 

The dot-level-depcndence of the current is presented 
in Fig. 2 for various bias voltages. In the linear response 
regime (V ~ 0), when the dot level passes through the 
Fermi energy of the electrodes, the resonant tunneling oc- 
curs and the current exhibits a typical Fano peak. When 
the bias voltage increases, the max-to-min distance be- 
comes larger. This can be understood as follows. The 
current peak is generated by the interference between 
electrons passing through two conduction paths and this 
takes place when the dot energy e g lies between the two 
chemical potentials. Since the chemical potential differ- 
ence is proportional to the bias voltage, the max-to-min 
distance of the current increases with the bias voltage. 
A similar deformation of the Fano peak as a function of 
the bias voltage was reported by Zhu et al. 13 

The AB oscillations at a fixed dot level are depicted 
in Fig. 3 for various bias voltages. The pattern of oscilla- 
tion changes when one of the chemical potentials passes 
through the dot level. In particular, an oscillation with 
half wavelength appears when \eV\ = 0.10 x ep with ep 
the Fermi energy at zero bias. 

As is known, 14 the two-terminal conductance G sat- 
isfies the Onsager relation G(ip) = G(—<p) and, thus, 
is an even function of <p. This implies that the conduc- 
tance takes its extremum at <p = or n. Ucda et al. 15 
pointed out that the conductance can take its maximum 
at ip ^ 0, 7r. It is interesting to note that these two cases 
are simultaneously observed in Fig. 3. When the bias volt- 
age increases, the maximum of the current at ip = splits 
into two maxima, which move to cp = ±ir in a continu- 
ous way. But the minimum of the current abruptly jumps 
from ip = 7r to <p = at the bias voltage where the oscil- 
lation with half wavelength appears. Similar AB oscilla- 
tion patterns can be seen in the linear response regime 
by varying the dot level. 16 

The strong bias-voltage-dependence of the AB oscil- 
lation suggests a possibility of a device sensitive to the 
magnetic field. Indeed, the current-bias-voltage charac- 
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Fig. 2. Nonlinear current at zero temperature as a function of 
the gate voltage applied to the QD for different bias voltages eV 
from below eV/e F =0.01, 0.02, 0.03, 0.06, 0.07, 0.085, 0.10, 0.15, 
0.25, 0.30. up is the Fermi energy at eV = 0. The AB phase is 
fixed to ip = and u 2 DW = 0.1. 
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Fig. 3. Phase dependence of the current. The gate voltage is fixed 
to the A-B line shown in Fig. 2. Other parameters are the same 
as those given in Fig. 2. 
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Fig. 4. Current versus bias voltage in units of e/h. The phase 
affects the current-bias- voltage characteristics: the phase (a) ip = 
(the solid line), (b) (p = n/2 (the dashed line) and (c) (p = n 
(the dotted line). 

teristics at a fixed dot level depend on <p (Fig. 4). 

4. Current Fluctuations 

Following Shimizu and Ueda, 17 we investigate the fluc- 
tuation of the charge q(r) = f Q T dt.J{t) leaving from the 
left electrode during the time interval r, where J(t) is 
the current operator from the left electrode defined just 
before eq.(4). Its fluctuation 

r^oo 7- 

can be evaluated with the aid of Wick's theorem and the 
two-point function (3) and is given by 

A/=(^) 2 £ rf e{T(e)X:F A (e)[l-F A (e)] 

+ T(e)[l-T(e)][F L (e)-F R (e)] 2 }. (7) 

The fluctuation AI consists of thermal and shot noises 
and, at zero temperature, only the shot noise survives. 

Fig. 5 shows a typical behavior of AI as a function of 
the dot level. One finds that AI has two dips, corre- 
sponding to the maximum and minimum of the current. 
Obviously, when the current is small, its fluctuation is 
small. And thus, the current minimum corresponds to the 
fluctuation dip. On the other hand, for one-dimensional 
perfect conductor, the current fluctuation is known to 
vanish 17 because all the available states are used for the 
transport and there remains no room for fluctuations. In 
the same way, when the current takes its maximum, the 
dot level is used for the transport and is not available for 
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fluctuations. Then, the fluctuation is suppressed and the 
dip appears. Similar results are reported in connection 
with the photo-assisted shot noise. 18 

When cp = 7r/2, the current is known to exhibit a sym- 
metric Breit-Wigner-type resonance peak as for the sys- 
tem without direct tunneling between electrodes (the sin- 
gle dot system) . This seems to indicate the absence of in- 
terference between electrons passing through two paths. 
However, it is not the case and the interference can be 
seen in the current fluctuation. Indeed, if the interfer- 
ence is absent, the fluctuation would be a sum of the 
fluctuation of the single dot system and constant fluctu- 
ation of the system without quantum dot. The curve (b) 
of the left figure of Fig. 6 depicts the fluctuation of the 
present system at ip = n/2, clearly which is not a sum of 
a constant function and the fluctuation of the single dot 
system (the curve (b) of the right figure of Fig. 6). The 
difference can be attributed to the interference between 
electrons passing through the two paths. 
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Fig. 5. The relationship between the current (a) and fluctuation 
(b) at low bias voltage range (eV/ep = 0.01), <p = 0, in units of 
(a) e/h and (b) (e/2nh) 2 . 
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Fig. 6. The current (a) and fluctuation (b) with the present sys- 
tem at ip = 7r/2 (left) and with the single dot system (right) in 
units of (a) e/h and (b) (e/2nh) 2 at eV/ep = 0.01 respectively. 

5. Conclusions 

In summary, exact NESS for the noninteracting model 
of an AB ring with a QD is explicitly constructed and 
bias-voltage-dependences of the Fano resonance and the 
AB oscillation are investigated. The AB oscillation pat- 
terns are shown to strongly depend on the bias voltage. 
This implies the nonlinear conduction sensitive to the 
magnetic field. The current fluctuations are studied as 
well and we show that, even when the interference can- 
not be seen in the average current, it contributes the 
current fluctuation. 
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Appendix: C*-algebraic approach to NESS 

In this approach, 10 one starts from a set of all (spa- 
tially) local observables, which form an algebra called the 
C*-algebra. Remind that, if the interaction range is finite, 
the commutator i[H, A] of the total Hamiltonian H and 
a local observable A is a local operator. Thus, both sides 
of the Heisenberg equation of motion hQ£ — i[H, A] are 
local and the time evolution A{t) can be treated within 
the framework of local observables. 

Physical states p are introduced by specifying aver- 
age values (A) p of all local observables A (i.e., the states 
are (positive) linear functionals over the algebra of lo- 
cal observables). Then, the long-term behavior of the 
system is naturally studied through the behaviors of 
the averages (A) p uy Then, for various systems, one can 
show the existence of the weak limits p± defined by 
lim.t^±oo(A) p / t j = (A) p± for any observable A. In this 
way, one can avoid the explicit use of infinite quantities 
such as the total Hamiltonian and total particle numbers, 
which are mathematically ill-defined, and one may real- 
ize the unidirectional state evolution in the weak sense. 

As mentioned in the text, our model is a typical ex- 
ample of the junction systems studied by Frolich et al. 9 
Since they rigorously show the existence of NESS and 
give its mathematical characterizations, we give a heuris- 
tic proof of the existence of NESS and its characteriza- 
tion presented in the text. 

Let Ho = Hl + Hr and define the M0ller operator 
fi = lim e iHt/h e -iH t/h = lim e -im/h e iH t/h^ thcn 

t — > — oo t — >+oc 

the incoming field operators are expressed as 

f3 kaX = lim e lHt / t k- tHot ^a k<T xe tH » t ^b- tHt / h = Qa klT xn 1 '. 

Suppose that the total Hamiltonian H docs not ad- 
mit bound states, then, since fi intertwines H with Hq: 
HVL = fiiJ and H has bound states represented by 
Ccr, fi should annihilate the states involving dot elec- 
trons and, thus, the product n = fit 17 is a projec- 
tion onto the states consisting only the electrons rep- 
resented by a ka \ (see the arguments in a standard 
text book on scattering theory 19 ). Therefore, one has 
n = (1 — cjcf)(l — ejej and it commutes with 
afeo-A. On the other hand, remind that the initial state 
can be written as pa oc p c cxp [— ^2\(H\ — fi\N\)/T\] 
where N\ is the total number operator of each 
electrode and p c is the initial dot state. Since po 



commutes with Hq, the steady state is given by 

p + = lim e -iHt/h e iH t/h poe -iH t/h e iHt/h_ Thig seemg 
t — >+oo 

to imply p + = fi/9 fit, but the normalization changes 
as in the wave function rcnormalization, and one has 
p+ = fi / 9ofi t /tr(fi / o fi t ) = ^/Oo^V( n )o, where (• • • ) 
stands for tr(- • • p ). Hence, for example, 

(^fel CTl Al Pk 2 <7 2 A 2 ^ fc 3 "3 As flk 4 <T 4 A 4 ) 

(7! X 1 Ua l 2 (7 2 A 2 ng fc3<T3 A3 ^4(74 A 4 g) 

<n> 

= ( a L<7 1 A 1 4 2 (7 2 A 2 a fc3(73A3afe4(74A4n)o/ (II)o 
= ( a l 1 (7 1 A i a l 2 (7 2 A 2 a fe3(73A3 a *:4(74A4)0 • 

The last equality follows from the fact that a k a\ and 
c a (hence, n) are statistically independent at the state 
po- Therefore, since the initial state po satisfies Wick's 
theorem for a k( r\, the same is true for p + and f3ka\ an d 
the two-point functions {Pl a \Pk'a'\') are given by (3). 
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